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Using Gibbs ensemble Monte Carlo simulations and density functional theory we investigate the 
fluid-fluid demixing transition in inhomogeneous colloid-polymer mixtures confined between two 
parallel plates with separation distances between one and ten colloid diameters covering the complete 
range from quasi two-dimensional to bulk-like behavior. We use the Asakura-Oosawa-Vrij model in 
which colloid-colloid and colloid-polymer interactions are hard-sphere like, whilst the pair potential 
between polymers vanishes. Two different types of confinement induced by a pair of parallel walls 
are considered, namely either through two hard walls or through two semi-permeable walls that 
repel colloids but allow polymers to freely penetrate. For hard (semi-permeable) walls we find that 
the capillary binodal is shifted towards higher (lower) polymer fugacities and lower (higher) colloid 
fugacities as compared to the bulk binodal; this implies capillary condensation (evaporation) of the 
colloidal liquid phase in the slit. A macroscopic treatment is provided by a novel symmetric Kelvin 
equation for general binary mixtures, based on the proximity in chemical potentials of statepoints 
at capillary coexistence and the reference bulk coexistence. Results for capillary binodals compare 
well with those obtained from the classic version of the Kelvin equation due to Evans and Marini 
Bettolo Marconi [J. Chem. Phys. 86, 7138 (1987)], and are quantitatively accurate away from 
the fluid-fluid critical point, even at small wall separations. However, the significant shift of the 
critical polymer fugacity towards higher values upon increasing confinement, as found in simulations, 
is not reproduced. For hard walls the density profiles of polymers and colloids inside the slit 
display oscillations due to packing effects for all statepoints. For semi-permeable walls either similar 
structuring or fiat profiles are found, depending on the statepoint considered. 

PACS numbers: 64.70.Ja, 82.70.Dd, 61.20.Gy 



I. INTRODUCTION 



Capillary effects that are induced by the confinement 
of a system are crucial to a variety of phenomena. An 
everyday example is the capillary rise, against gravity, 
of the meniscus of a free gas-liquid interface at the wall 
of a container that encompasses the fluid. The contact 
angle at which the gas-liquid interface hits the wall is 
described by Young's equation 1], 7igcos0 = 7„g — 7^1, 
where the relevant interfacial tensions are those between 
the coexisting liquid and gas phase, 7ig, the wall and the 
gas phase, 7wg, and the wall and the liquid phase, 7wi. 
When the contact angle is finite, < 9 < n, the liquid 
phase partially wets the wall; for 6 = a, macroscopic 
layer of liquid grows between the wall and the gas phase 
far away from the wall, hence the liquid completely wets 
the wall; for — tt a correspondingly reversed situation 
occurs: a macroscopic layer of gas grows between the 
wall and the liquid phase far away from the wall, hence 
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the wall is completely "wet" by the gas phase and drying 
occurs. All these phenomena are driven by the influence 
of a single wall on the fluid. 

Different, but related effects occur under confinement 
of fiuids in pores, e.g. between two parallel planar walls. 
The phase diagram of such a confined system can dif- 
fer significantly from that in bulk 0, |^ . Depending on 
the nature of the interactions between fiuid particles and 
the walls, the bulk gas with chemical potential /i < /isat, 
where /isat is the value at saturation, can condense inside 
the pore and form a dense confined liquid phase. This 
capillary phase transition is accompanied by a jump in 
the adsorption isotherm at a given value of fi < /isat- 
Confinement may lead to stabilization of a phase that 
is unstable (or at least metastable) in bulk for a given 
statepoint. As opposed to capillary condensation the op- 
posite effect, referred to as capillary evaporation, is also 
feasible: A fluid with chemical potential fj, > /Xsat forms a 
liquid in bulk, but evaporates inside the capillary. A sim- 
ple, yet powerful, way to quantitatively describe capillary 
phase transitions is based on the Kelvin equation, derived 
from a macroscopic treatment of the thermodynamics of 
the inhomogeneous system. For capillaries with planar 
slit-like geometries [j, it predicts that liquid-gas coexis- 
tence inside the slit will occur at /x = /^sat + 7ig cos{9)/h, 
where h is the separation distance between the two walls. 
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Hence, contact angles Q < 6 < tt /2 correspond to cap- 
illary condensation, while 7r/2 < 6* < tt corresponds to 
capillary evaporation. 

In this work we investigate capillary phenomena oc- 
curring on a mesoscopic scale using a simple model for 
a mixture of sterically-stabilized colloidal particles and 
non-adsorbing polymer coils confined between two paral- 
lel planar walls. Using colloids as model systems offers 
advantages over molecular substances due to easy exper- 
imental access to the relevant length and time scales, and 
the possibility of using real space techniques like confo- 
cal microscopy. Our work is intended to stimulate such 
investigations in order to gain further insights into capil- 
lary phenomena. The topology of the bulk phase diagram 
of colloid-polymer mixtures depends on the polymer size 
(as compared to the size of the colloids) and the polymer 
concentration. At sufhciently high polymer concentra- 
tion and for sufficiently large polymer, the bulk system 
demixes into a colloidal liquid phase that is rich in colloids 
(and dilute in polymers), and a colloidal gas phase that 
is dilute in colloids (and rich in polymers) 0, 0, B ■ 
This phenomenology is similar to gas-liquid coexistence 
in simple fluids with the polymer fugacity Zp playing 
the role of inverse temperature, and hence governing the 
strength of interparticle attractions. However the phase 
separation in colloid-polymer mixtures is purely entropy- 
driven, due to an effective depletion interaction between 
the colloidal particles, that is mediated by the polymer 
[Tol ITU . A similar mechanism induces an effective deple- 
tion interaction between a colloidal particle and a pla- 
nar hard wall 12i|. For a review on the properties of 
colloid-polymer mixtures in contact with a hard wall, 
we refer the reader to Brader et al. \13t . Recently the 
wall-fluid interfacial tensions 0, lisl Il6| and the contact 
angle of the free interface and a hard wall ^3 j the 
gas-liquid interfacial tension 0, 0| have been studied 
theoretically and with computer simulations. Moreover 
complete wetting of the colloidal liquid phase at a hard 
wall is found, close to the critical point, in both simula- 
tion [ISIHllll and theory [lllljl. Complete wetting 
of the colloidal liquid is also observed in experimental 
realizations of colloid-polymer mixtures in contact with 
glass substrates j23i IHl or glass substrates that possess 
the same coating with an organophilic group as the col- 
loidal particles do [2^. On the other hand, experiments 
(27| with polymer-grafted substrates (of the same chem- 
ical nature as the dissolved polymers) showed that the 
contact angle is larger than 7r/2. Although the structure 
of the polymer coating was not studied in detail, Wijting 
et al. VZq expect the polymers to form a fluffy layer with 
a distribution of loops and tails. Such a layer is known to 
diminish the depletion interaction between colloidal par- 
ticles and the substrate 28] . Wessels et al. |15j modeled 
this type of substrate using a semi-permeable wall, that 
is completely penetrable to the polymers but act like a 
hard wall on the colloids. Using DFT they predict com- 
plete drying. More recently the behavior of the mixtures 
in contact with porous walls was investigated, and both 



wetting of the surface and drying into the porous matrix 
depen ding on the precise path in the phase diagram were 
found [23 ■ Porous walls are wet by the liquid phase, 
with a transition from partial to complete wetting at a 
polymer fugacity Zp almost independent on the porosity 
of the wall. 

Although much research has been devoted to the un- 
derstanding of the behavior of the mixture in contact 
with a single wall, few studies have been reported on 
the effect of confinement on the phase behavior of the 
mixture. Computer simulations and theory were used 
to study porous matrices. When impenetrable to both 
components they were found to induce capillary conden- 
sation, whereas matrices penetrable to the polymers, but 
not to the colloids, induce capillary evaporation [30ll3l|. 
Furthermore laser induced confinement |3^ was found 
to induce capillary condensation. The only experimental 
result on capillary condensation that we are aware of is 
that by Aarts et al. 25] . These authors found capillary 
condensation for colloid-polymer mixtures confined in a 
wedge formed by a glass bead placed on a glass substrate. 

The aim of the present paper is to determine the phase 
behavior and the structure of colloid-polymer mixtures 
confined between two parallel planar walls for a complete 
range of wall separation distances. We focus on the fluid 
part of the phase diagram, and present results for hard 
walls and semi-permeable walls, obtained from computer 
simulations and density functional theory. A preliminary 
account of this work was given in Refs. |33l l34j . Here we 
supply all relevant details of the simulation method. We 
also investigate the structure of the mixtures inside the 
pore, and derive and test a generalized Kelvin equation 
for binary mixtures confined in slit-like pores and com- 
pare its predictions with our simulation results. 

II. MODEL 

Consider a binary mixture of Nc hard spheres of di- 
ameter (Jc representing the colloidal particles and iVp 
non-interacting spheres of diameter tjp representing the 
(ideal) polymers in a volume V at temperature T. This 
is a simple model for a colloid-polymer mixture as the 
interaction between sterically-stabilized colloids can be 
made such that it resembles closely that of hard spheres, 
and a dilute solution of polymers in a 0-solvent is weakly 
interacting. Colloids and polymers interact via a hard- 
sphere-like potential, as the polymers are excluded from 
a centre-of-mass distance (ctc -I- o'p)/2 from the colloid, 
where Up = 2Rg, and Rg is the radius of gyration of the 
polymer coils. We treat the solvent as an inert contin- 
uum. This so-called Asakura-Oosawa-Vrij (AOV) model 
[IE HE US is defined by the pair potentials 

. . J 00 if Rij < CTc 
Vcc(nij) - <^ Q otherwise, ^ ' 

where Rij = \Ri — Rj\ is the distance between two col- 
loidal particles, with Ri the position of the centre-of-mass 
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a) b) 




FIG. 1: (a) Illustration of the AOV model of hard sphere col- 
loids (gray circles) of diameter CTc and ideal polymers (trans- 
parent circles) of diameter cTp confined between parallel hard 
walls of area A and separation distance H. Colloids behave as 
hard spheres, polymers cannot penetrate colloids, and poly- 
mers may freely overlap. The walls are impenetrable to both 
components. The 2-axis is perpendicular to both walls, and 
the origin is located in the middle of the slit, (b) Same as 
in (a) but for walls that are penetrable for polymers and im- 
penetrable for colloids (semi-permeable walls) of area A and 
separation distance He- Polymers are confined between paral- 
lel hard walls (continuous line) with separation distance Hp. 
This is a model for substrates at distance Hp coated with a 
polymer brush (not shown) of thickness {Hp — He) /2; the free 
distance between both polymer brushes is He- Solute poly- 
mers (open spheres) are able to penetrate the brush (dashed 
lines) but not the substrate; the brush acts like a hard wall 
to the colloids. 

of colloid i, 

vpp{nj)=0, (2) 

where Vij = \ fi—fj \ is the distance between two polymers, 
with fi the position of the centre-of-mass of polymer i and 

^ ■' [0 otherwise, 

where \Ri —rj \ is the distance between colloid i and poly- 
mer J. 

The size ratio q — (Jp/ac is a geometric parameter that 
controls the range of the effective depletion interaction 
between the colloids. We denote the packing fraction 
by Vk = (7rCT^A^fe)/(6F), with fc = c, p for colloids and 
polymers, respectively. As alternatives to rjp, we use as a 
thermodynamic variable the colloidal fugacity Zc, or the 
polymer reservoir packing fraction r]^ that satisfies the 
(ideal gas) relation 

Vl = Jf^p^p- (4) 



The hard walls are modeled such that neither polymers 
nor colloids can penetrate the walls. The wall-particle 
potential acting on species fc = c, p is 

fO \f-{H-ak)/2<z<{H-ak)/2 
Vwk{z) - I ^ otherwise, 

(5) 

where z is the coordinate normal to the walls, and H is 
the wall separation distance. We define the volume of 
the system as = AH, where A is the (lateral) area of 
the confining walls. Fig.QJa) shows an illustration of the 
model. 

Semi-permeable walls are defined by the external po- 
tential 

, , r if -{Hk - ak)/2 < z < [Hk - afe)/2 
«wfel^j - I oo otherwise, 

(6) 

where He is the wall separation distance that the colloids 
experience, while Hp is the wall separation distance that 
the polymers feel; the latter can be interpreted as the 
distance between substrate walls. We define the volume 
as = AHc- Fig. QJb) shows an illustration of this 
model. Any choice of Hp > Hc + 2ap leads to decoupling 
of the (inner) colloid and (outer) polymer wall, since the 
ideal gas of polymers cannot mediate correlations from 
the substrate to the interior of the system. 

III. SIMULATION METHOD 

We perform Gibbs Ensemble Monte Carlo (GEMC) 
simulations to determine the phase behavior of the AOV 
model in bulk and confined between two parallel planar 
walls. The number of colloids Nc, the number of poly- 
mers A'^p, and the volume V are kept fixed and are divided 
into two separate subsystems a and b with volume and 
V^, respectively, with the constraints that V = + 
, Nc = + and Np = + N^. The two subsys- 
tems are allowed to exchange both particles and volume 
in order to satisfy the conditions for phase equilibrium, 
i.e., equality in both phases of the chemical potentials of 
the two species and of the pressure. The method can be 
applied to bulk as well as to confined systems [ssL Is^ , 
and will be described briefly below; for more details, we 
refer the reader to Ref. [4(1 |. 

In more detail, the bulk equilibrium conditions be- 
tween two coexisting phases require equal temperature 
T, equal chemical potential /i; for each species i (z = c, p 
in our case) and equal pressure P. In our model, the 
temperature T is irrelevant; because of the hard-core na- 
ture of the interaction potentials it acts only as a scaling 
factor setting the thermal energy scale. In the GEMC 
method, the two coexisting phases are simulated simul- 
taneously in two separate (cubic) boxes with standard 
periodic boundary conditions. The acceptance proba- 
bility for a trial move to displace a randomly selected 
particle is 

r = min{l, exp[-/3(C/„ew - C/oid)]}, (7) 
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where 



exp[-/3([/„ew - C^old)] = 



1 if /?[/, 
if f3U, 



new 
new — 



(8) 



with C/now and [/qm the energy of the new and old gen- 
erated configuration, and /3 = 1/(/cbT) the inverse tem- 
perature with ks being the Boltzmann constant. Note 
that Uoid = as configurations with UoM = oo are ex- 
cluded by a vanishing Boltzmann factor (regardless of 
the temperature). Transfer of (single) particles between 
the two boxes is used to satisfy equal chemical potential 
for both species. We select at random which subsystem 
is the donor and which is the recipient as well as the 
species (colloid or polymer) of the particle transfer. Sub- 
sequently, a specific particle is randomly selected in the 
donor box and transferred to a random position in the 
recipient box with probability 



V = min < 1 , 



aT/b 



■exp[-/3(f7„ew-t/oid)] L (9) 



where i = c, p is the species of the selected particle, 
is the volume of the donor box and is the volume of 
the recipient box. 

In addition, volume changes of the boxes are used to 
satisfy the condition of equal pressure in both subsys- 
tems. We hence perform a random walk in ln(T4/Vfc) 
with the acceptance probability 



P = min{l,7eexp[-/3(t/„ 



where 



7^ 



old 



b 

new 



^old 



Uou)]}: 



N°+N° + l 



(10) 



(11) 



and with the condition that the total volume V 



old 



''old 



l^icw is constant , where the sub- 



script old (new) marks quantities before (after) the trial 
displacement. 

Phase equilibria in confined geometries, like e.g. inside 
slit pores as considered here, can be determined by either 
determining the adsorption isotherm in simulation or by 
employing the GEMC simulation method extended to a 
slit pore [39]. In the first method, the pore is coupled 
to a reservoir of bulk fiuid. The adsorption of the fiuid 
inside the pore is then measured at constant temperature 
for different bulk densities of the reservoir, i.e. the chem- 
ical potentials of the various species are fixed. A jump 
in the adsorption isotherm corresponds to capillary con- 
densation or evaporation in the slit pore. However, this 
method is inaccurate for determining phase coexistence 
due to hysteresis. 

In order to determine the binodal lines, we hence em- 
ploy the GEMC method adapted to a slit pore [s^ : Two 
separate simulation boxes are simulated, one containing 
the confined gas and one containing the confined liq- 
uid. Each box has periodic boundary conditions in both 



directions parallel to the walls. Capillary coexistence 
implies equality of temperature, chemical potentials for 
both species as well as equality of the wall-fluid interfacial 
tension. Fulfilling the first two conditions is performed 
similar to the procedure for determining bulk phase equi- 
libria, i.e. using particle displacements with acceptance 
probability given by Eq. Q and particle transfers with 
acceptance probability given by Eq. We satisfy the 
third requirement of equal wall-fiuid interfacial tension in 
both phases by exchanging wall area (and hence volume) 
between the two boxes, while fixing both the wall sepa- 
ration distance H in each box as well as the total lateral 
area of both boxes + = A is constant, with A'^ 
and A^ the area of the subsystems a and b, respectively. 
A random walk in \n{A^/A^) is then performed with an 
acceptance probability given by 



r = min{l,7eexp[-/3(t/„ 



C/old)]} 



(12) 



with 



n^i^^] . (13) 



A 



old 



^old 



We determine the fugacities of the colloids and of the 
polymers by applying the GEMC version of the particle 
insertion method l4ll 



ya,b 



a.b 



1 



■exp[-/3A[/] 



c, p 



(14) 



where a, b label the two simulation boxes and AU is the 
energy defined by the acceptance rule of Eq. ((HJ. We de- 
termine the densities of the two coexisting phases by sam- 
pling histograms of the probability density P(?7c,?7p) to 
observe the packing fractions rip and ijc for the polymers 
and colloids, respectively. The two maxima of P{ric,rip) 
correspond to the coexisting packing fractions in the ther- 
modynamic limit |4l| . Statistical uncertainties of the 
sampled quantities were determined by performing three 
or four independent sets of simulations. The standard de- 
viation of the results from the simulation sets was used 
as the error estimate. 



IV. KELVIN EQUATION FOR BINARY 
MIXTURES 



In this Section, we derive expressions for the shift in 
chemical potentials for the gas-liquid binodal of a binary 
mixture confined between parallel plates assuming knowl- 
edge of bulk quantities like the bulk coexisting densi- 
ties, the chemical potentials at bulk coexistence, and the 
liquid- wall and gas- wall interfacial tensions, 7wi and 7wg, 
respectively. Using a macroscopic picture we employ the 
grand potential of the mixture confined in a slit of two 
parallel walls with area A and separation distance h be- 
tween the two walls 52j 

il(^c, Mp) = Ahuj[iic, Mp) + 2A7(/^c, Mp, h), (15) 
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FIG. 2: Phase diagram of the AOV model with size ratio 
q = cTp/cTc = 1 as a function of the coUoidal chemical potential 
fic and the polymer chemical potential fip. The bulk binodal 
from free volume theory (thick continuous curve) is shown 
along with isobaric lines (thin continuous curves) for two bulk 
reference points. The thick dashed curves indicate the normal 
paths for the same reference points. The thin dashed curves 
are an illustration of the possible prediction of the Kelvin 
equation for the capillary binodals. For clarity we omitted 
the paths with constant polymer chemical potential. 



where /ic and fip are the chemical potentials of colloids 
and polymers, respectively, to is the bulk grand potential 
density (per unit volume), and 27 is the interfacial ten- 
sion of the mixture and the plates at distance h. For large 
/i, we can approximate the latter quantity by twice the 
interfacial tension of the mixture in contact with a single 
wall, e.g., 27wQ(/ic, Mp): where a=l,g are for the liquid 
and the gas phase, respectively. The aim is to predict 
the grand potential in the capillary given the knowledge 
of the thermodynamics of the bulk at coexistence, i.e. 
at a statepoint specified through the bulk values of the 
chemical potentials /i* and /i*. Hence, we can reexpress 
the chemical potentials of both species for the confined 
fluid as /Xc = /i* + A^c and /ip = //* + A/ip. Quanti- 
ties at coexistence carry a superscript a = 1, g, where 1, g 
indicate liquid and gas respectively. 

The thermodynamic relations for the bulk densities of 
the colloids and polymers, pc and pp, respectively, read 



mers, Tc and Fp, are given by 



F. = 



(17) 



Using Eqs. (|16|l and IjlTfl . we can perform a Taylor ex- 
pansion of the right hand side of Eq. H15|l 



Ah 



^^(McMp) 



2 r 



7wa+CA/ic-t-F^A/J,(18) 



where the bulk densities, p" and and the (excess) 
adsorptions, F" and Fp, are evaluated in one of the co- 
existing phases a = 1, g at the statepoint given by /x* and 

We next consider the capillary to be at phase coexis- 
tence, i.e. we might envisage two capillaries in contact 
with each other, one being filled with gas, the other be- 
ing filled with liquid. Phase equilibrium between both 
capillaries implies equality of the chemical potential of 
each species and of the grand potential in both phases. 



(19) 



where ftg and fii is the grand potential of the gas and 
the liquid phase, respectively. Using the approximation 
(El in (PJ) yields 



(F^-Ff) A^i, 



Pp-p|,-TK-r|,)\A^^p 



(7wi(A*c'Aip) -7wg(Aic'Mp))- (20) 



In the limit h is large compared with the microscopic 
lengths and that the adsorptions remain finite, we can 
neglect the terms proportional to {T\ ~ Tf)/h. Hence, 
Eq. Ipljl simplifies to 



(Pc - P?)Aa*c + (Pp - P^)A/ip = ^(^ 



7wi — 7wj 



(21) 



which constitutes a single equation for the two unknown 
shifts in the chemical potentials, A/ic and A/ip. In or- 
der to obtain a closed system of equations one requires 
another assumption. We will present three different ap- 
proaches in the following. Each approach can be viewed 
as a different choice of an "optimal" bulk reference state 
as illustrated in Fig. |21 A comparison with the numerical 
results will be presented in Sec. IV El 



duj 
dpc 



Pp 



duj 
dpp 



(16) 



while the (excess) adsorption of the colloids and the poly- 



A. Constant Polymer Chemical Potential 

The bulk reference state is chosen to be at the same 
polymer chemical potential. This requirement is similar 
to imposing the same temperature in a simple substance. 
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The optimal bulk reference state with the same chemical 
potential of polymers leads to 

A^p = 0. (22) 

Eqs. (|2H) and H22|l readily yield a one-component Kelvin 
equation 



The derivation presented here yields the same results as 



A 2 "^"1 

A^c = T — 



7« 



(23) 



Clearly this result can be obtained with less labour by di- 
rectly dealing with the effective one-component system of 
colloids interacting with a depletion potential, and rather 
serves as a illustration for the validity of the reasoning 
leading to Eq. (|^ . 



B. Constant Pressure 

This reference state was used by Evans and Marini 
Bettolo Marconi in Ref. Q. The bulk reference state is 
chosen to possess the same pressure. In order to derive 
a corresponding condition consider the (finite difference 
version of the) Gibbs-Duhem relation 

{S/V)AT-AP + p,A^I, + ppA^ip = 0, (24) 

where S is the entropy. Clearly, for athermal systems 
AT = 0. We use (|24|l to compare the statepoint at capil- 
lary coexistence with the reference statepoint at bulk co- 
existence that possesses the same pressure, i.e. AP — 0, 
hence 



A/ip 



;A/ic 



a 



,1. 



(25) 



As shown in Fig. [5]the constant pressure paths (for which 
the Gibbs-Duhem relation is a good approximation close 
to the bulk binodal) have a discontinuity at bulk coexis- 
tence. To predict the capillary phase behavior for state- 
points that lie on the gas side of the bulk gas-liquid bin- 
odal (hence to predict capillary condensation), the cor- 
rect reference state is the coexisting gas phase, and we 
use (|25|l with a = g. Together with 1)21(1 we obtain the 
classic result for capillary condensation, 



A^c 
A^p 



^ 7wi - 7wg 
2 7wi - 7wg 

^Pp - {pI/p^)pV 



(26) 



Alternatively, predicting the capillary phase behavior for 
statcpoints that lie on the liquid side of the bulk gas- 
liquid binodal (and to predict capillary evaporation), we 
choose the coexisting liquid phase as a reference. Eq. H25|l 
with a = 1 together with Eq. H21|l leads to the following 
result for capillary evaporation 



Apc 
A/ip 



7wi - 7w 



hpl- [pl/p\,)Pv ' 
2 7wi - 7wg 

hpl- {pI/pDpI' 



(27) 



npr 

given in Ref. where a single capillary is considered 
in contact with a bulk reservoir. We note that this pro- 
cedure leads to two different equations for the two phe- 
nomena of capillary condensation and evaporation. 



C. Normal Path Relation 

As a third and novel alternative we choose the reference 
state such that the statepoint of interest lies on a path in 
the phase diagram that crosses the bulk liquid-gas bin- 
odal in perpendicular (or normal) direction in the plane 
of chemical potentials. This implies optimizing proximity 
in the pc — Pp plane between both statepoints. 

The corresponding relation can readily be established 
starting from the relation of the slope of the bulk binodal 
to the difference in coexisting densities. 



dpc 
dpp 



Pc-Pc 



(28) 



Hence a path normal to the bulk binodal is given by 



dpc 
dpp 



Pc - Pc 



imal 



Pi. 



pp 



from which we deduce the finite-difference version 



A/ip 



Pp -Pj 



Pc - Pc 

which, together with Eq. (|21|) yields 



(29) 



(30) 



A/ic 
App 



2. 

^(7wi 
^(7wi 



7wg) 



7wg) 



ipl-p'c)' + {p'p-p'p)' 

Pp - Pp 

ipl-p'c)' + ip'p-p'p)' 



(31) 



valid both for capillary condensation and evaporation. 
As shown in Fig. |21the normal paths are symmetric with 
respect to the gas and liquid side of the bulk binodal. 
The somewhat subtle choice of the reference state, which 
is different for the two phenomena in Sec. IV b, resulting 
in using either Eq. H26|) or Eq (|27|l is now avoided. The 
procedure described here leads to one equation, which 
can be used for both phenomena. This might be advan- 
tageous for mixtures where the corresponding phases are 
less obvious, like e.g. confined liquid crystals. 



V. RESULTS 

In this section, we focus on colloid-polymer mixtures 
with a size ratio g = 1, a previously well-studied and 
also experimentally accessible case. We have performed 
GEMC simulations with 4 x 10* steps discarding the ini- 
tial 10* steps for equilibration. The acceptance proba- 
bility of particle displacement was kept around 10% to 
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FIG. 3: Bulk phase diagram of the AOV model with size ratio 
q = (Jp/(Jc = 1 as a function of the colloidal packing fraction 7)c 
and the polymer reservoir packing fraction rf^ . Coexistence is 
along the horizontal tie lines (not shown). Shown are GEMC 
simulation results for volumes V — SOOOcrj? (circles) and V = 
1250ct? (diamonds), as well as the simulation results from 
Ref. pCJ] (squares). We also display the result of the fit of 
Eqs. (1321 and 13311 (solid line) to the simulation data to and 
the binodal obtained from DFT, or, equivalently, free volume 
theory (dashed line). 



20%, the acceptance probability for volume exchanges 
was around 50%, while the acceptance probability for 
transfer of particles was strongly dependent on the state- 
point and varied between 50% and 5% for the polymers 
and from 10% to less than 0.1% for the colloids. The 
probability of a successful particle transfer typically de- 
creases with increasing particle density. The statistical 
error of the fugacities serves as a good indicator for the 
efficiency of the particle transfer between the two boxes; 
an unreasonably high error indicates that the simulation 
results are unreliable due to too low particle transfer 
probabilities. Using this indicator the maximum pack- 
ing fraction that could be reached in our simulations was 
?7c ~ 0.34. 

Fig. |51 shows the bulk phase diagram obtained from 
simulations with a simulation box of volume V = SOOOcr^. 
Simulation runs with smaller system sizes display neg- 
ligible finite size effects. In the rjc-rf^ representation, 
the shape of the binodal is similar to that of a simple 
fluid upon identifying rf^ with inverse temperature. The 
tielines, connecting the coexisting phases, are horizontal 
(not shown) as rf^ possesses the same value in the two co- 
existing phases. We have checked that our results agree 
well with those obtained by Dijkstra and van Roij 
who performed simulations of an effective one component 



FIG. 4: Capillary phase diagram of the AOV model with size 
ratio q — ap/ac — 1 confined between parallel hard walls with 
separation distance H/(Jc=2, 3, 5, 10 and cx) as a function 
of colloid packing fraction rj^ and polymer reservoir packing 
fraction rjp. Shown are results from simulation (symbols) and 
DFT (continuous lines). 

system, which was obtained by formally integrating out 
the degrees of freedom of the polymers in the partition 
function. 

To estimate the location of the critical point, we fitted 
the binodals using the scaling law 

and the law of rectilinear diameter 

where TyJ, is the colloid packing fraction of the coexist- 
ing liquid phase, ryf is the colloid packing fraction of the 
coexisting gas phase and the subscript crit indicates the 
value at the critical point. A and B are two free pa- 
rameters determined from the fit, and (3 = 0.32 is the 
three-dimensional Ising critical exponent. We used the 
standard functional form of the two laws ^3 but replaced 
the temperature by the inverse of the polymer reservoir 
packing fraction. The continuous curve in Fig.|31is the re- 
sult of the fit of Eqs. and The fit is remarkably 
good, but we point out that this only gives an estimate 
of the critical point. To get a more precise value of the 
critical packing fractions, it would be necessary to carry 
out simulations in a region of the phase diagram much 
closer to the critical point than is possible with GEMC. 
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(77p)crit 


(T7c)crit 


DFT inlUt 


DFT (r7c)crit 


oo 


0.86(1) 


0.117(2) 


0.638 


0.103 


10 


1.00(1) 


0.124(1) 


0.670 


0.120 


5 


1.09(1) 


0.123(2) 


0.710 


0.111 


3 


1.27(2) 


0.116(3) 


0.815 


0.100 


2 


1.76(2) 


0.119(1) 


1.044 


0.091 



2.5 



TABLE L Capillary critical points of the AOV model be- 
tween parallel hard walls with separation distance H/ac — 
oo, 10, 5, 3, and 2 obtained from the fit of Eqs. 1)32^ and II33|I . 
and from DFT. 



In Fig. O we also compare our results with those ob- 
tained from density functional theory (DFT). We use the 
approximation for the Helmholtz excess free energy for 
the AOV model as given in ^j. For given external po- 
tential, the density functional is numerically minimized 
using a standard iteration procedure. The discrepan- 
cies between theory and simulation can be understood by 
considering that the DFT for homogeneous (bulk) fluid 
states of the AOV model is equivalent to the free vol- 
ume theory of Lekkerkerker et al. 0|. Dijkstra et al. 
showed that this theory is equivalent to a first order 
Taylor expansion of the free energy around — 

/3F(7Ve, V, 77P = /3F(7Ve, V, = 0) + 

^''^^..rv^5/3^(^c,v^,(r/;,)0\ 

' (%)' 



[ d{ril)'[- 



6 



it < Vfrcc > ,;,=o , (34) 



^3 ''P 



neglecting terms 0{{'q^p)'^) and where {Vircc)-qi=o is the 
free volume available for the polymer in the pure hard 
sphere reference system. It is evident in Fig. O that the 
theory (dashed curve) performs better at high rjc where 
the system is so crowded that it resembles the reference 
hard sphere system and {Vbee) (rii)' =q - {Vbee) (r,i)' =rii- 
For very small rjc the free volume is close to the total 
volume of the system V for both the pure hard sphere 
reference system and the actual mixture. For high rf^, the 
gas-liquid coexistence is very broad and quantitatively 
well-predicted by the theory. The critical point of the 
AOV mixture for size ratio g = 1 is in the region of 
rjc ~ 0.1 and the theory underestimates the critical value 
of T^p. Furthermore, the discrepancy in location of the 
critical point arises from the mean field critical exponent 
of the theory against the 3D Ising critical exponent of 
the simulation 19]. 



A. Phase diagrams 

First we present the results for the colloid-polymer 
mixtures confined between two smooth, planar hard walls 
at distance H . In Fig. 21 we show a set of phase diagrams 
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FIG. 5: The same as Fig. 2] but as a function of colloid fugac- 
ity Zcot and polymer reservoir packing fraction rf^. Statepoint 
A is in the gas region of the phase diagram for wall separa- 
tion distances H/ac ~ 00, 10, and 5 and in the liquid region 
of the phase diagram for wall separation distances H/ac = 3 
and 2. The arrow indicates the direction of the binodal shift 
upon increasing confinement (decreasing values of H/ac) of 
the mixture between parallel hard walls. 



for H/ac — 00 (bulk), 10, 5, 3, and 2, in the rjc-rjp repre- 
sentation. Upon decreasing the plate separation distance 
H/ac, the critical value of rf^ shifts to higher values, in 
accordance with the decrease in critical temperature of 
simple fluids. The theoretical binodals agree well with 
those from simulation, except close to the critical point. 
The theory underestimates for all plate separations the 
critical value of 77p, as it does for the bulk system 
We observe that the deviation increases upon decreasing 
H/ac- In Tab. n we show the critical packing fractions 
obtained from the fit of Eqs. (|32|l and (|33|l and from 
DFT. We used the 3-dimensional Ising critical exponent 
for all wall separations. Although recent studies [ijl sug- 
gests a critical behavior for small wall separations that is 
neither three-dimensional, nor two-dimensional, the dif- 
ference is likely to be negligible at the level of precision 
of our GEMC simulations. 

In Fig.|Slwe show a set of phase diagrams for H/ac = 
00, 10, 5, 3, and 2 in the Zc-rf^ representation. The co- 
existence gap in colloid packing fractions collapses to a 
line since two phases at coexistence possess the same col- 
loid fugacity. Note that the system is in the gas phase 
for fugacity Zc < z^^*, while it is in the liquid phase for 
Zc > Zc*^', where z^*** denotes the colloid fugacity at bulk 
coexistence. Statepoint A is gas-like for H/ac > 5, but 
is liquid- like for H/ac < 3. Hence, planar slits with 
H/ac < 3 are filled with liquid phase, while the bulk 
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(77p)crit 


(T7c)crit 


DFT inlUt 


DFT (r7c)crit 


oo 


0.86(1) 


0.117(2) 


0.638 


0.103 


10 


1.09(2) 


0.13(1) 


0.660 


0.075 


4 


1.11(4) 


0.10(3) 


0.699 


0.076 


2 


1.29(1) 


0.11(4) 


0.818 


0.092 



TABLE II: Capillary critical points of the AOV model be- 
tween parallel semi-permeable plates with separation distance 
H/ac = oo, 10, 4, and 2 obtained from the fit of Eqs. (I32II and 
CT . and from DFT. 



reservoir is in the gas phase, proving the occurrence of 
capiUary condensation upon reducing the width of the 
sht that is in contact with a bulk gas. 

We now turn our attention to coUoid-polymer mixtures 
confined between two smooth, planar semi-permeable 
walls at distance He- In Fig.Elwe show a set of phase dia- 
grams for Hc/cTc — oo, 10, 4, and 2 in the 7?c-?7p represen- 
tation. Upon increasing the confinement (via reduction 
of Hc/(Tc), the critical value of rj^ shifts to higher values. 
The trend is similar to the behavior of the slit with hard 
walls, although we find a smaller shift of the critical point 
(see Tab. EJ. In Fig.|7|we show a set of phase diagrams 
for H/fJc = 00,10,4, and 2 in the Zc-rj^ representation. 
Statepoint B is liquid-like for ffc/cc > 10, but is gas-like 
for Hc/(Jc < 2. Hence, planar slits with Hc/(Tc < 2 are 
filled with gas, while the bulk reservoir is in the liquid 
phase, indicating the occurrence of capillary evaporation. 




FIG. 6: Capillary phase diagram showing the gas-liquid bin- 
odal for the AOV model with q — a-^ja^ — 1 between 
two parallel semi-permeable walls with separation distances 
i?c/o"c=2, 4, 10, oo as a function of colloid packing fraction 
77c and polymer reservoir packing fraction i]p. Shown are re- 
sults from simulations (symbols) and DFT (lines). 



B. Structure at coexistence 

We next analyze the density profiles of both species 
at capillary coexistence of gas and liquid phases. Such 
fluid states are translationally invariant against lateral 
displacements, and the density distributions (of both 
species) depend solely on the (perpendicular) distance 
from the walls. We compare theoretical and simulation 
results for coexistence states at the same values for rjp. 
In practice, we have used the result for rj^ from the simu- 
lations, and have calculated the corresponding DFT pro- 
files. The value for Zc used in the DFT calculations was 
adjusted according to the respective theoretical capillary 
binodal. Recall that the quantitative differences in re- 
sults for the capillary binodals from simulation and the- 
ory are small. 

The results for a slit of hard walls are shown in Fig. |S1 
The colloidal profiles in the liquid phase display strong 
layering at either wall. For wall separation H / cjc = 10, 
at 77;; = 1.23 ± 0.01, and z^a^ = 66.3 ± 0.1 these oscil- 
lations decay to flat, bulk-like behavior in the center of 
the slit. For smaller wall separations, namely H/dc — 5, 
at Tf^ = 1.39 ± 0.01, and z^al = 99.4 ± 0.3, H/a^ = 3, 
at if^ = 1.68 ± 0.01, and z^al = 128 ± 2, and H/a^ ^ 2, 
at ril = 2.23 ± 0.01, and z^al = 228 ± 10, we observe 
the presence of 5, 3, and 2 well-defined layers of parti- 
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FIG. 7: The same as Fig. |S| but as a function of colloid fu- 
gacity Zca^ and polymer reservoir packing fraction rf^. The 
arrow indicates the direction of the shift of the binodal upon 
increasing confinement of the mixture (decreasing values of 
Hc/(Jc) between two parallel semi-permeable walls. 
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FIG. 8: Density profiles of the coexisting liquid (left) and 
gas (right) phases of the AOV model confined between 
parallel hard walls with separation distance: a) H/(Jc=10, 
at 77^=1.23(l),zccr?=66.3(l), b) H/crc=5 



=99.4(3), c) H/(7e=3, 



at r;J,=1.68(l), 

.3_ 



at 

.3_ 



?7p=1.39(l), 
2ccrj = 128(2), and 

d) H/crc=2, at »7;;=2.23(1), 2ccrj^=228(10). Shown are re- 
sults from simulations for the density profiles of the colloids 
(crosses) and polymers (circles), along with results from DFT 
(full lines). 



cles, respectively. Although their density is much lower, 
the polymers in the liquid phase display similar behavior. 
The layering is weaker, but we can observe that a maxi- 
mum in the colloid profile corresponds also to a maximum 
in the polymer profile. This result suggests that for such 
low concentrations polymers behave as hard spheres as 



packing effects are concerned. In the gas phase, for wall 
separations H/ac = 10 and 5, we find strong adsorption 
of the colloids at both walls, and a tendency of the poly- 
mers to desorb from the walls. In the center of the slit 
almost no colloids are present and the polymers display 
flat density profiles with a packing fraction very similar to 
the polymer reservoir packing fraction. For wall separa- 
tion distances H/ac = 3 and 2, we observe an almost flat 
polymer density profile, while the density of the colloids 
is very low throughout the slit. Different from the liquid 
profiles, a maximum in the colloidal profile corresponds 
to a minimum of the polymer profiles. 

Fig. El displays density profiles for the slit of semi- 
permeable walls. For wall separation distance of H^/crc = 
10, at = 1.39±0.01, and z^al = 167±5, we clearly ob- 
serve for the liquid state points the presence of a gas layer 
between the wall and the liquid phase centered in the slit. 
For wall separation Hc/<Jc = 4, at if^ = 1.35 ± 0.01, and 
ZcCrl — 318±10 the gas layers at the walls disappear and 
indications of layering effects appear. For wall separation 
i/c/(Tc = 2, at 77^; = 1.61±0.01, and z^al = 2950±200 the 
colloid density profile displays very significant peaks at 
both walls. Moreover, we do find layering at larger wall 
separations, for statepoints well inside the liquid phase. 
We will discuss this in more detail in the next section. In 
the gas phase, the density of colloids is very low through- 
out the slit, while the polymer density profile is almost 
flat with a packing fraction close to the polymer reservoir 
packing fraction. 

The comparison between DFT and simulation indi- 
cates good agreement of results from both approaches. 
Differences in structure can be traced back to differences 
in the phase diagrams. For fixed rf^ the DFT predicts 
higher colloid densities and smaller polymer densities, 
and these differences are reflected by the discrepancies 
in the profiles. For wall separations where the statepoint 
is close to the critical point the agreement is worse, es- 
pecially close to the walls. Such discrepancies between 
simulation and DFT were previously reported in study- 
ing the wall-fluid tension and the adsorption of colloid- 
polymer mixtures [l^ . 



C. Structure ofF-coexistence 

We next consider density profiles for a fixed state- 
point off-coexistence. For slits with hard walls we chose 
statepoint A (77p=1.49 and ln(zcO-j?)=4.6) of Fig.|Sl that 
lies in the stable gas region of the bulk phase diagram. 
We carried out simulations for wall separation distances 
H/oc — 10, 5, 3, and 2. Fig. 1101 shows that for wall sepa- 
rations H/ac — 10 and 5 the slit is filled with gas. How- 
ever, for wall separation oi H/ac = 3 we observe that 
the capillary fills with liquid. Hence, for this particular 
statepoint, the critical wall separation distance for capil- 
lary condensation lies between 3 and 5 colloid diameters, 
consistently with the findings of Sec. IV Al Reducing the 
wall separation to H/ac = 2, the liquid phase remains 
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uid fills slits with wall separations Hc/uc = 10, and 4, 
while for Hc/cTc = 2 the slit is filled with gas. This is an 
indication of capillary evaporation consistent with the 
findings of Sec. IV Al The liquid phase is characterized 
by structureless polymer profiles, and a layering of col- 
loids for both Hc/<Jc — 10, and 4. Note that we did not 
observe such layering for the statepoints at coexistence 
(see previous section). Simulation and theory are in good 
agreement at all wall separations. 
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FIG. 9: Density profiles of the coexisting liquid (left) and gas 
(right) phases of the AOV model confined between two par- 
allel semi-permeable walls with varying separation distance; 
a) Hc/crc = 10, at 77;;=1.39(1), 2ca?=167(5), b) Hc/crc=4, at 
ryp = 1.35(l), 2ccr,^=318(10), and c) Hc/crc=2, at 77;; = 1.61(1), 
Zccrc=2950(200). Shown are results from simulations for the 
density profiles of the colloids (crosses) and polymers (circles), 
and results from DFT (full lines). 



FIG. 10: Density profiles of the AOV model with q — 1 be- 
tween parallel hard walls with varying separation distance: 
a) H/ac^lO, b) Hc/ac^ 5, c) H^/a^=2,, and d) Hc/(J,=2 at 
statepoint A of the phase diagram of Fig. |21 i.e., for poly- 
mer reservoir packing fraction r^p = 1.49 and colloid fugacity 
ln(2ccrc) = 4.6. Shown are results from simulations for the 
density profiles of the colloids (crosses) and polymers (cir- 
cles), and results from DFT (full lines). 



stable. The density profiles in the gas phase possess ad- 
sorption peaks in the colloid profile, and corresponding 
desorption peaks in the polymer profiles. In the liquid 
phase we observe strong layering of the colloids, and, to 
a lesser extent, of the polymers. The agreement between 
simulation and DFT results is good. The differences seem 
to be related to the vicinity of statepoint A to the critical 
point. The critical point for wall separation H/gc = 3, 
and 5, is much closer to statepoint A than the critical 
points for wall separations Hj = 10, and 2, where we 
find a better agreement between simulation and theory 
profiles. 

We next discuss the structure of the mixture inside the 
slit with semi-permeable walls (Fig. Ill|l . We fix the fu- 
gacities of both species to those at statepoint B (77p=1.60 
and ln(0cfTj?)=7.3), see Fig. [71 The statepoint B is in the 
liquid part of the bulk phase diagram phase. The liq- 



D. Two dimensional limit 

We next analyze the dimensional crossover from three 
to two spatial dimensions by reducing the distance of 
the hard walls towards 7J/ctc — > 1. The two-dimensional 
system encountered for Hja^ = 1 is identical to a two- 
dimensional mixture of colloidal hard disc and ideal poly- 
mer discs. Two-dimensional mixtures were previously 
studied with both theory '45', '4^ , and experiments |47j |. 
For H very close to CTp the polymer reservoir packing 
fraction scale as 
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FIG. 11: Density profiles of the AOV model between two par- 
allel semi-permeable walls with varying separation distance: 
a) Hc/ffc^W, b) //c/o"c=4, and c) Hc/ac=2 at statepoint B of 
the phase diagram of Fig. |7| i.e., for polymer reservoir pack- 
ing fraction rip — 1.60 and colloid fugacity ln(zcO"c) — 7.3. 
Shown are results from simulations for the density profiles of 
the colloids (crosses) and polymers (circles), and results from 
DFT (fuU lines). 



We eliminate the divergence using scaled variables for 
the polymer reservoir packing fraction rf^{H — a-p) / H and 
for the colloidal fugacity Zc{H — ac)/H. Effectively, we 
map the three-dimensional system with packing fractions 
rii = TTcrfNij i&AH) to the two-dimensional system with 
packing fractions rji = ■KafNi/{4:A), where z = c, p. 

In Fig. I12r a) we plot the phase diagrams in the scaled 
Vp ~ Vc representation and we observe that the binodals 
for H/ac — 1-01 and H/ac = 1.005 are superimposed, 
demonstrating that this is a reliable estimate for the 
binodals of the 2-dimensional system. The compari- 
son with a two-dimensional DFT kg equivalent to a 
two-dimensional free volume theory |45j indicates poorer 
agreement as in the three-dimensional case. The discrep- 
ancy in the critical polymer reservoir packing fraction is 
very substantial. We also find that in this representa- 
tion the binodals of the three-dimensional system of the 
slits collapse over the bulk binodal, indicating a scaling 
of the critical value of r]p as (?7p)crit ~ ^/{H ~ cTp). In 
Fig. ITW b) we see that the collapse of the binodals onto a 
master curve in the scaled rj^— scaled Zc representation 
is not as good as in the other representation, moreover 
the sequence of binodals is inverted, indicating that the 
critical value of the colloid fugacity scales differently. 




(H-a,)/H) 



FIG. 12: Phase diagrams showing the binodals for the AOV 
model between parallel hard walls with separation distance 
H/ac = cxD, 10, 5, 3, 2, 1.01, and 1.005. Shown are results from 
simulation (symbols), from three-dimensional DFT (dashed 
curves), and from DFT for the AOV model in two dimen- 
sions (full lines), a) The gas- liquid binodal as a function of 
the scaled variable rjp{H — ap)/H and r^c; b) The gas- liquid 
binodals as a function of the scaled variables Zcffc (H — ac)/H 



and ril{H^Gp)/H. 



E. Kelvin equation 

In this section, we compare the simulation results with 
the predictions of the Kelvin equations that we derived in 
section Hvl First, we address the relationship of the pa- 
rameter h, we used in the Kelvin equations, to our model 
parameter H . There are two, a priori equivalent, choices 
that we investigate, namely h = H, and h = H — cjc- 
Second, the Kelvin equations need, as a input, the dif- 
ference of the gas and liquid tensions at the wall inter- 
face. Since this data are not readily available we will 
assume, in the case of capillary condensation, the rela- 
tion — = Jig , strictly valid only in the complete 
wetting regime, to hold at all state points considered. 
Likewise, for the capillary evaporation case, we assume 
7u)i ~ Iwg — lig , valid in the complete drying regime to 
hold at all state points. For the liqu id- g as interfacial 7ig 
tension we use DFT data from Ref. 

Figs. Il^f a) and ll^f c) display the simulation results 
for the hard wall slit together with the predictions of the 
Kelvin equations H22|) , H2t)|l , and (|31|l for h — H — (Jc, and 
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h = H respectively. The Kelvin equation (|22|l . derived 
from the path with constant polymer chemical poten- 
tial, is superimposed at all separation distances with the 
Kelvin equation l|26|) derived using the constant pressure 
path. This is consistent with the observation of Aarts 
and Lekkerkerker [2^. Now we can offer an alternative 
explanation. As shown in Fig [3 for the capillary con- 
densation case the bulk in the gas phase (polymer rich 
phase) , and the path with constant polymer chemical po- 
tential is almost equivalent to the constant pressure path. 
In other words, the bulk reference point for Eq. H22I) and 
Eq. (|26|l is very similar in the case of capillary conden- 
sation. The Kelvin equation 1)31(1 derived using a normal 
path, predict a smaller shift with respect to Eq. ((26)1 . To 
estimate the error introduced by the complete wetting 
approximation {'^wg — Iwi = 7/g)j we show few points 
(filled diamonds) predicted by Eq. I(26() using the actual 
difference in wall tensions, as published in our previous 
work 33] . 

The Figs. IT^ b) and IT^ d) display the simulation re- 
sults for the semi-permeable wall slit together with the 
predictions of the Kelvin equations 1(22(1 . 1(271) . and l(31() 
for h = H — Gc, and h = H respectively. The prediction 
of Eq. l|77)l and (|!TT|l are superimposed at all separation 
distances considered. This is surprising since they are 
derived from very different "paths", as shown in Fig. |21 

We can conclude that the Kelvin equation we derived 
using a novel approach, gives predictions that are con- 
sistent, and in quantitatively agreement with the predic- 
tion of the classic equation. In addition our equation 
is the same for capillary condensation and evaporation, 
and the choice of reference state, which is different for 
both phenomena, is avoided. As the Kelvin equation is 
based on macroscopic arguments, is surprising that we 
find reasonable quantitative agreement for nearly two- 
dimensional systems. The shift of the critical polymer 
fugacity towards higher values upon increasing confine- 
ment, as found in simulations, is not reproduced, because 
the Kelvin equation is entirely based on properties of the 
(semi) infinite system. Finally, the two choices of h we 
presented give essentially the same results for wall sepa- 
rations H/dc as small as 4. 



VI. CONCLUSIONS 

In conclusion, we have studied the effect of strong con- 
finement provided by two parallel walls on the phase be- 
havior and structure of model mixtures of colloids and 
polymers of size ratio q = \. The densities of the gas 
and liquid phases at coexistence, as well as the chemi- 
cal potentials were computed by GEMC simulations and 
DFT. Two different models of confining walls were in- 
vestigated: i) Hard walls, impenetrable to both colloids 
and polymers, as a model for glass walls in contact with 
colloid-polymer mixtures were found to stabilize the liq- 
uid phase for statepoints that lie in the gas part of the 
bulk phase diagram; this effect is referred to as capillary 
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FIG. 13: Phase diagram of colloid-polymer mixtures in the 
colloidal chemical potential log(zcO-j?) and the polymer reser- 
voir packing fraction rf^ representation, a) Results for slits 
with hard walls, and h = H — Oc- b) Results for the slit with 
semi-permeable walls , and h = H — ac- c) Results for the 
slits with hard walls, and h = H . d) Results for the slit with 
semi-permeable walls, and h — H. In a) and c) the simula- 
tion results (symbols) are compared with the prediction of the 
Kelvin equations 1221 , 12611 , and 13111 assuming complete wet- 



ting {-yu 



The filled diamonds are results from 



Eq. 12611 using the difference in wall tensions 7„g — 'y^i taken 
from Ref. 33]. In b) and d) the simulation results (symbols) 
are compared with the prediction of the Kelvin equations 1221 , 
1271 . and 1311 assuming complete drying (7^,; — -y^g = 719) • 
For the liquid-gas interfacial tension 7ig we used DFT data 
from Ref. ,48]. Few lines are superimposed, see text for the 
explanation. 



14 



condensation, ii) Semi-permeable walls, impenetrable to 
colloids but penetrable to polymers, could be experimen- 
tally realized using polymer coated substrates [23| ■ If the 
coating density is not too high the polymer brushes can 
act as impenetrable for colloids while being penetrable for 
polymers. We find that the effect of semi-permeable walls 
is to stabilize the gas phase for statepoints in the liquid 
part of the bulk phase diagram; this effect is referred to 
as capillary evaporation. Both capillary evaporation and 
condensation are consistently predicted by GEMC sim- 
ulations and DFT. The differences between simulations 
and DFT in the vicinity of the critical point are con- 
firmed for bulk mixtures and were found to be larger for 
the confined mixtures. The differences reach a maximum 
in the limit of two-dimensional colloid-polymer mixtures. 

We have studied the structure of the mixture between 
parallel walls by measuring the density profiles in the 
direction normal to the confining walls. For the liquid 
phase, rich in colloids and poor in polymers, we found 
layering of colloids with an oscillation period roughly 
equal to the diameter of the particles for all wall sep- 
arations and statepoints considered. In the case of semi- 
permeable walls, the structure, i.e. the layering of colloids 
and the adsorption or desorption of gas layers at the semi- 
permeable walls, depends strongly on the statepoint and 
on the lengthscale of the confinement. For the gas phase, 
rich in polymers and poor in colloids, we found flat poly- 
mer profiles with moderate desorption of polymers from 
the hard walls. We found that density oscillations for 
colloids and polymers are correlated in the liquid phase 
and anti-correlated in the gas phase. This can be un- 
derstood by the following argument. In the liquid phase 
the fraction of polymers is small and a polymer is always 
surrounded by other colloids to which it interacts with an 
hard-core potential. Clearly the polymer structure must 
be similar to that of the colloids. In the gas phase the 



fraction of polymers is large with respect to the colloids 
and the entropy is increased by segregation of colloids, 
so if a region is locally denser in polymer than the av- 
erage polymer density, it will be more dilute in colloids. 
The comparison between simulation and DFT is overall 
good, with small differences in the vicinity of the critical 
point. Our findings of capillary condensation for hard 
walls and capillary evaporation for semi-permeable walls 
are consistent with the experimental findings of Aarts 
et al. nil and Wijting et al. jHH^. Nevertheless the 
comparison is only qualitative. Experiments with better 
controlled geometries are needed to relate directly to our 
predictions. The use of the surface force apparatus, for 
example, should be possible for colloid-polymer mixtures. 
Simulations and theory should proceed to more realistic 
models for colloid-polymer mixtures [i^ Isol l5lj |. For 
example, the inclusion of excluded volume interactions 
between polymer coils has given accurate results for the 
bulk phase diagram _4QJ . The use of realistic models for 
confined mixtures should help understanding the surface 
phase behavior of real colloid-polymer mixtures. Work 
along this line is in progress. 
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